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These lectures provide an introduction to lattice methods for nonperturbativc stud- 
ies of Quantum Chromodynamics. Lecture 1 (Ch. 2) is a very vanilla introduction 
to lattice QCD. Lecture 2 (Ch. 3) describes examples of recent lattice calculations 
relevant to fixing the parameters of the CKM matrix. 



1 Introduction 

The latticeB regularization of QCD has been a fruitful source of quahtative and 
quantitative information about QCD for many years, especially when combined 
with Monte Carlo simulation. Lattice methods are presently the only way we 
know how to compute masses and matrix elements in the strong interactions 
beginning with the Lagrangian of QCD. My goal in these lectures is to give 
enough of an overview of the subject that you will be able to make an intelligent 
appraisal of a lattice calculation. 

The first lecture will describe how to put QCD on a lattice. This is a long 
story with a lot of parts. Lattice QCD is full of technicalities, but I will try 
to make the discussion physical. In Lecture Two I will discuss lattice calcula- 
tions of matrix elements which are needed to convert experimental numbers to 
predictions for the CKM matrix. These calculations have a lot of ingredients: 
a typical one starts with a particular choice of discretization and simulation 
algorithm, and a choice of operators whose matrix elements are appropriate 
for one's measurement. After the lattice number is computed, it may have to 
be converted into a number in a continuum regularization scheme (like MS), 
which will involve some kind of perturbative or nonperturbative matching cal- 
culation. Finally, it might have to be extrapolated in quark mass, to some 
physical quark mass value or the the chiral limit. None of these parts are sim- 
ple or obvious (or more precisely, most of the time the simple and obvious idea 
doesn't work very well). Hopefully you will find the physics in the calculations 
more interesting than the tables of numbers which result. 



2 Basics of Lattice QCD 
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2. 1 Lattice Variables and Actions 



All quantum field theories must be regulated in order to control their ultraviolet 

divergences while calculations arc performed. The lattice is a space-time cutoff 
which eliminates all degrees of freedom from distances shorter than the lattice 
spacing a. As with any regulator, it must be removed after renormalization. 
Contact with experiment only exists in the continuum limit, when the lattice 
spacing is taken to zero. The lattice is a unique regulator compared to the 
ones you might already know. Other regularization schemes are tied closely to 
perturbative expansions: one calculates a process to some order in a coupling 
constant; divergences are removed order by order in perturbation theory. The 
lattice, however, is a nonperturbative cutoff. Before a calculation begins, all 
wavelengths less than a lattice spacing arc removed. 

All regulators have a price. On the lattice we sacrifice all continuous 
space-time symmetries but preserve all internal symmetries, including local 
gauge invariancc. This preservation is important for nonperturbative physics. 
For example, gauge invariance is a property of the continuum theory which is 
nonperturbative, so maintaining it as we pass to the lattice means that all of 
its consequences (inchiding current conservation and renormalizability) will be 
preserved. The bill is paid when we take the lattice spacing to zero and try to 
recover what we have left out. 

Let's begin by thinking about a lattice version of scalar field theory. One 
just replaces the space-time coordinate a;^ by a set of integers {x^ = an^, 
where a is the lattice spacing). Field variables (^(a;) are defined on sites (/'(a;„) = 
(j)n, The action, an integral over the Lagrangian, is replaced by a sum over sites 

S= f d'^xL^a'^Yl^iM- (1) 

n 

and the generating functional for Euclidean Green's functions is replaced by 
an ordinary integral over the lattice fields 

Z= [([[dcp^)e^. (2) 

Gauge fields are a little more complicated. They carry a space-time index ^ in 

addition to an internal symmetry index a (^J^(a;)) and arc associated with a 
path in space a;^(s): a particle traversing a contour in space picks up a phase 
factor 

ip P(exp ig / dXf,A^)i; = U{s)ip{x). (3) 

J s 
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P is a path-ordering factor analogous to the time-ordering operator in ordinary 
quantum mechanics. Under a gauge transformation g, U{s) is rotated at each 
end: 

Uis) ^ g-\x^{s))U{s)g{x^{0)). (4) 

These considerations led Wilson i to formulate gauge fields on a space-time 
lattice, in terms of a set of fundamental variables which are elements of the 
gauge group G living on the links of a four-dimensional lattice, connecting 
neighboring sites x and x + afi: U^{x), with C/^(a; + /i)^ = Ufj,{x) 

U^,{n)^exp{igaT-A;in)) (5) 

for SU{N). {g is the coupling, the vector potential, and is a group 
generator) . 

Under a gauge transformation link variables transform as 

U^ix)^V{x)U^ix)V{x + ii)^ (6) 

and site variables as 

^{x) V{x)ip{x) (7) 

so the only gauge invariant operators we can use as order parameters are matter 
fields connected by oriented "strings" of U's 

i!{xi)U^(xi)Uf,{xi+ fi)...ij(x2) (8) 

or closed oriented loops of U's 

Tr . . . U^ix)U^{x + A) . . . ^ Tr . . . U^{x)V^ix + fi)Vix + (i)U^,{x + /*).... (9) 

An action is specified by recalling that the classical Yang-Mills action in- 
volves the curl of F^u. Thus a lattice action ought to involve a product 
of J7^'s around some closed contour. Gauge invariance will automatically be 
satisfied for actions built of powers of traces of U's around arbitrary closed 
loops, with arbitrary coupling constants. If we assume that the gauge fields 
are smooth, we can expand the link variables in a power series in gaA'^s. For 
almost any closed loop, the leading term in the expansion will be proportional 
to F'^^. This is not a bug, it is a feature. All lattice actions are just bare 
actions characterized by many bare parameters (coefficients of loops). In the 
continuum (scaling) limit all these actions are in the same universality class, 
which is (presumably) the same universality class as QCD with any regulariza- 
tion scheme, and there will be cutoff-independent predictions from any lattice 
actions which are simply predictions of QCD. 
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Let's hold that thought while we do an example: 

The simplest contour has a perimeter of four links. In SU{N) 

^=^EEl^eTr (l-C/4n)f/.(n + A)C/t(„ + ^)f/t(„)). (iq) 

n ^>v 

This action is called the "plaquette action" or the "Wilson action" after its 
inventor, is the bare lattice coupling, whose associated cutoff is a. The 
lattice parameter (3 = 2N/g^ is often written instead of g"^ = 47ras. 

Let us see how this action reduces to the standard continuum action. Spe- 
cializing to the U(l) gauge group, and slightly redefining the coupling, 

5 = 4 ^ ^ Re {l-exp{iga[Ai_,{n)+A^{n+fi)-Ai,{x+i>)-A^{n)])). (11) 

^ n > 

The naive continuum limit is taken by assuming that the lattice spacing a is 
small, and Taylor expanding 

A^{n + i>) = Ai,{n) + ad^Ai,{n) + . . . (12) 

so the action becomes 

/3S=^^^l-Re (exp(i5a[a(a,A^ - d^A,) + ©(a^)])) (13) 

= i«'EE^^^.^^ + -- (14) 

= \j d^^Fl (15) 
transforming the sum on sites back to an integral. 

2.2 Numerical Simulations 

In a lattice calculation, like any other calculation in quantum field theory, 
we compute an expectation value of any observable F as an average over a 
ensemble of field configurations: 

{V) = ^ jmeM-SW)- (16) 

We do this by Monte Carlo simulation: we construct an ensemble of states 
(collection of field variables), where the probability of finding a particular con- 
figuration in the ensemble is given by Boltzmann weighting (i. e. proportional 
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to exp(— 5*). Then the expectation value of any observable F is given simply 
by an average over the ensemble: 



1 ^ 

<r>^f ^-^r[0W]. (17) 



i=l 



As the number of measurements N becomes large the quantity T will become 
a Gaussian distribution about a mean value, our desired expectation value. 
The idea of essentially all simulation algorithms B is to construct a new config- 
uration of field variables from an old one. One begins with some initial field 
configuration and monitors observables while the algorithm steps along. After 
some number of steps, the value of observables will appear to become indepen- 
dent of the starting configuration. At that point the system is said to be "in 



equilibrium" and Eq. 17 can be used to make measureraents. 

Dynamical fermions are a complication for QCD The fermion path 
integral is not a number and a computer can't simulate fermions directly. 
However, one can formally integrate out the fermion fields. For nf degenerate 
fermion flavors 

„ "/ 

Z= [dU][diP][d^p] exp(-/3S'G(t/) - V^M(C/)V') (18) 

i=l 

= j[dU]idetM{U))"f exp{-pSiU)). (19) 
The determinant introduces a nonlocal interaction among the C/'s: 

Z = j [dU] exp{-/3S{U) - n/Trln(M(C/))). (20) 

Generating configurations of the f/'s involves computing how the action changes 
when the set of C/'s are varied. Typically, this involves inverting the fermion 
matrix M{U) {d\og M / dM — M~^). This is the major computational prob- 
lem dynamical fermion simulations face. M has eigenvalues with a very large 
range- from 2t: down to mqa- and in the physically interesting limit of small 
niq the matrix becomes ill-conditioned. At present it is necessary to compute 
at unphysically heavy values of the quark mass and to extrapolate to = 0. 
(The standard inversioiLtechnique today is one of the variants of the conju- 
gate gradient algorithm B. ) This tremendous expense is responsible for one 
of the "standard" lattice approximations, the "quenched" approximation. In 
this approximation the back-reaction of the fermions on the gauge fields is 
neglected, by setting nf = in Eq. |l^. Valence quarks, or quarks which 



5 



appear in observables, are kept, but no sea quarks. No one knows how good 
an approximation this is, in principle. In practice it works very well for spec- 
troscopy. The only way we know how to test it is to compare simulations in 
the quenched approximation with those from full QCD. 

2.S Spectroscopy Calculations 

Masses are computed in lattice simulations from the asymptotic behavior of 
Euclidean-time correlation functions. A typical (diagonal) correlator can be 
written as 

C{t) = {0\O{t)O{0)\0). (21) 

Making the replacement 

0{t) = e^*Oe-^* (22) 
and inserting a complete set of energy eigenstates, Eq. |2l| becomes 

Cit)=J2\{0\O\n)\^e'^-K (23) 

At large separation the correlation function is approximately 

C(t) ~ |(0|O|l)|2e-^i* (24) 

where Ei is the energy of the lightest state which the operator O can create 
from the vacuum. Fig. ^ shows an example of this. If the operator does not 
couple to the vacuum, then in the limit of large t one hopes to to find the mass 
El by measuring the leading exponential falloff of the correlation function. If 
the operator O has poor overlap with the lightest state, a reliable value for the 
mass can be extracted only at a large time t. In some cases that state is the 
vacuum itself, in which Ei — 0. Then one looks for the next higher state-a 
signal which disappears into the constant background. This is hard to do. 

Most of the observables we are interested in will involve valence fermions. 
Let's suppose we wanted to measure the mass of a meson. Then we might take 

X 

where 

J{x,t) =i>{x,t)T^p{x,t) (26) 

and r is a Dirac matrix. The intermediate states |n) which saturate C{x,t) 
are the hadrons which the current J can create from the vacuum: the pion. 
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Figure 1: An obviously very nice looking lattice correlator and its fit. Periodic boundary 
conditions convert the exponential decay into a hyperbolic cosine. 

for a pscudoscalar current, the rho, for a vector current, and so on. Now we 
write out the correlator in terms of fermion fields 

C{t) = ^(0|^i(a;,t)«ry^,(a;,i)«V^fe(0,0)''r«VKO,0)^|0) (27) 

X 

with a Roman index for spin and a Greek index for color. We contract creation 
and annihilation operators into quark propagators 

{Omj{x,trMO,Ofm = G^^{x,t;0,0) (28) 

so 

C{t) = ^ TrG(a;, t; 0, 0)rG(0, 0; x, t)T (29) 

X 

where the trace runs over spin and color indices. Baryons are constructed 
similarly. A good way to think about these correlators is by using a sort of 
Feynman-diagram language which keeps track of the valence quark lines but 
ignores all the gluons and sea quarks. 

2.4 The Continuum Limit 

When we define a theory on a lattice the lattice spacing a is an ultraviolet cutoff 
and all the coupling constants in the action are the bare couplings defined with 
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respect it. When we take a to zero we must also specify how the coupUngs 
behave. The proper continuum hmit comes when we take a to zero holding 
physical quantities fixed, not when we take a to zero fixing the couplings. 

On the lattice, if all quark masses are set to zero, the only dimensionful 
parameter is the lattice spacing, so all masses scale like 1/a. Said differently, a 
lattice calculation produces the dimensionless combination am{a). One can de- 
termine the lattice spacing by fixing one mass from experiment. Then all other 
dimensionful quantities can be predicted. Imagine computing some masses at 
several values of the lattice spacing. (Pick several values of the bare parame- 
ters and calculate masses for each set of couplings.) Our calculated mass ratios 
will depend on the lattice cutoff. If the lattice spacing is small enough, the 
typical behavior will look like 

(ami(a))/(aTO2(a)) = mi(0)/rn2(0) + 0{mia) + 0{{miaf) + . . . (30) 

The leading term does not depend on the value of the UV cutoff, while the 
other terms do. The goal of a lattice calculation is to discover the value of 
some physical observable as the UV cutoff is taken to be very large, so the 
physics is in the first term. Everything else is an artifact of the calculation. 
We say that a calculation "scales" if the a— dependent terms in Eq. ^ are zero 
or small enough that one can extrapolate to a = 0, and generically refer to all 
the a— dependent terms as "scale violations." 

We can imagine expressing each dimensionless combination am{a) as some 
function of the bare coupling(s) {g{a)}, am = f{{g{a)}). As a ^ we must 
tune the set of couphngs {g{a)} so 

lim —f{{g{a)}) constant. (31) 

a-^Q a 

From the point of view of the lattice theory, we must tune {g} so that correla- 
tion lengths 1/ma diverge. This will occur only at the locations of second (or 
higher) order phase transitions. In QCD the fixed point is gc = so we must 
tune the coupling to vanish as a goes to zero. 

One needs to set the scale by taking one experimental number as input. A 
complication that you may not have thought of is that the theory we simulate 
on the computer is different from the real world. For example, the quenched 
approximation, or for that matter QCD with two flavors of degenerate quarks, 
almost certainly does not have the same spectrum as QCD with six flavors 
of dynamical quarks with their appropriate masses. Using one mass to set 
the scale from one of these approximations to the real world might not give a 
prediction for another mass which agrees with experiment. 
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(The glass is always half empty.. .Ip the strong coupling limit, lattice reg- 
ularized. QCD automatically confines □ and chiral symmetry is spontaneously 
broken So unless there is some kind of phase transition as the bare cou- 
plings are tuned to take the cutoff away, which probably doesn't happen, we 
are working with a confining theory without doing anything special.) 

Today's QCD simulations range from 16^ x 32 to 32'^ x 100 points and 
run from hundreds (quenched) to thousands (full QCD) of hours on the fastest 
supercomputers in the world. The cost of a Monte Carlo simulation in a box 
of physical size L with lattice spacing a and quark mass ruq scales roughly as 

{-)H-r\—?-' (32) 

a a niq 

where the 4 is just the number of sites, the 1-2 is the cost of "critical slowing 
down" -the extent to which successive configurations are correlated, and the 
2-3 is the cost of inverting the fermion propagator, plus critical slowing down 
from the nearly massless pions. Thus it is worthwhile to think about how 
to do the discretization, to maximize the value of the lattice spacing. The 
thing to keep in mind is that the lattice action is just a bare action defined 
with a cutoff. No lattice discretization is any better or worse (in principle) 
than any other. Any bare action which is in the same universality class as 
QCD will produce universal numbers in the scaling limit. However, by clever 
engineering, it might be possible to devise actions whose scaling behavior is 
better, and which can be used at bigger lattice spacing. 

An example Q of a test of scale violations is shown in Fig. |[ The x axis is 
the lattice spacing, in units of a quantity ri , which is defined through the heavy 
quark potential: rfdV{r)/dr\r-^ — 1.0, about 0.4 fm. The plotting symbols are 
for different kinds of discretizations. The flatter the curve, the smaller the 
scale violations. 

The simplest organizing principle for "improvement" is to use the canon- 
ical dimensionality of operators as a guide. Consider the gauge action as an 
example. If we perform a naive Taylor expansion of a lattice operator like the 
plaquette, we find that it can be written as 

1 - ^Re TrUpiaq = roTrF^^ + a^[n E^. "^^D^F^.D^F^, + 

r2Y.,u.'^^D^F,,D^F,„ + 
r3E^..T^^-Df.F^aD,F,,] + 

+0(a4) (33) 

The expansion coefficients have a power series expansion in the coupling, rj — 
Aj + g^Bj + . . . and the expectation value of any operator T computed using 
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Figure 2; Lattice calculations of the (a) rho and (b) nuclcon mass, interpolated to the point 
m^ri = 0.778, as a function of lattice spacing. 



the plaquette action will have an expansion 

(T(a)) = (r(0)) + 0(a) + Oig'a) + ... (34) 

Other loops have a similar expansion, with different coefficients. Now the idea 
is to take the lattice action to be a minimal subset of loops and systematically 
remove the a" terms for physical observables order by order in n by taking 
the right linear combination of loops in the action, S = ^jCjOj with c 
c'j + g^Cj + . . .. This method was developed by Symanzik and co-workers 
in the mid-80 's. 

Ordinary perturbation theory (expansions in the bare lattice coupling 
are not very convergent, but cleverjirescriptions for definitions of couplings E. 
or nonperturbative tuning method£j have been quite successful in developing 
improved lattice actions. 

2.5 Relativistic Fermions on the Lattice 

Defining fermions on the lattice involves yet another problem: doubling. Let's 
illustrate this with free field theory. The continuum free action is 



S = / d^x[iPixh^d^ijix) + m^ix)i;ix)]. (35) 



10 



One obtains the so-called naive lattice formulation by replacing the derivatives 
by symmetric differences: 

Sr''" = E V5«^(V'n+M - ^n-^) + V^nV-n. (36) 



The propagator is: 



G{p) = smp^a + ma) = — — — — (37) 

2^ sm p^a + m^a^ 



We identify the physical spectrum through the poles in the propagator, at 
Po = iE: 

sinh^ Ea — sin^ pja + m^a^ (38) 
i 

The lowest energy solutions are the expected ones at p = (0,0,0), E ~ ±m, 
but there are other degenerate ones, at p = (tt, 0, 0), (0, tt, 0, ), ... (tt, tt, tt). As 
a goes to zero, the lightest excitations of the spectrum, the ones whose energy 
is 0(1), not 0(l/a), are the relevant ones, and there are sixteen of these, in all 
the corners of the Bruilloin zone. Thus our action is a model for sixteen light 
fermions, not one. This is the famous "doubling problem." j— . 

In fact, associated with the "doubling problem" is the Nielsen- Ninomayallj 
theorem, which says that no lattice action can be undoubled, chiral, and have 
couplings which extend over a finite number of lattice spacings (ultralocality) . 
However, there are three ways to get two out of three. They are 
(a) Wilson Fermions (undoubled, nonchiral, ultralocal) 

We can alter the dispersion relation so that it has only one low energy 
solution. The other solutions are forced to i? ~ 1/a and become very heavy 
as a is taken to zero. The simplest version of this solution, called a Wilson 
fermion, adds an irrelevant operator, a second-derivative-like term 

= -^Y. - SV-n + V^n-M) - arV^i^V (39) 

to 5'"°'"'^ . The parameter r = 1 is almost always used and is implied when one 
speaks of using "Wilson fermions." 

There are two dimension-five operators which can be added to a fermion 
action. The Wilson term is just one of them. The other dimension-five term 
is a magnetic moment term 

Ssw - ^-^^j{x)af_,uFfj,uij{x) (40) 
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Figure 3: The "clover term" 



and if both terms are included, their coefficients can be tuned so that there are 

SO(ag^) lattice artifacts. This action is called the "Sheikholeslami-Wohlert" 
or "clover" action because the lattice version of is the sum of paths 
shown in Fig. ^. 

Wilson-type fermions contain an explicit chiral-symmetry breaking term. 
This causes a lot of bad things to happen. The most obvious is that the 
zero bare quark mass limit is not respected by interactions; the quark mass 
is additively renormalized. The value of bare quark mass niq which the pion 
mass vanishes, is not known a priori before beginning a simulation; it must be 
computed. This is done in a simulation involving Wilson fermions by varying 
iriq and watching the pion mass extrapolate quadratically to zero as ~ 
•niq — niq. It actually turns out that this is a worse problem than you would 
think: the Dirac operator _D on a gauge configuration could develop a real 
eigenmode A at minus the bare quark mass you dialed into the program. Then 
D + m would be non-invertible! Other nasty things happen (operator mixing, 
ge the next section) and people argue about how serious they are in practice 

(b) Staggered or Kogut-Susskind Fermions (chiral, doubled, ultralocal) 

In this formulation one reduces the number of fermion flavors by using 
one component "staggered" fermion flelds rather than four component Dirac 
spinors. The Dirac spinors are constructed by combining staggered fields on 
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Figure 4: An example of flavor symmetry breaking in an improved staggered artion. The 
different 7's are a code for the vario»Si pseudoscalar states. Data are from Ref. O. For an 
explanation of the splitting, see Ref. tj. 



different lattice sites. Staggered fermions preserve an explicit chiral symmetry 
as VTLq — > even for finite lattice spacing, as long as all four flavors are degen- 
erate, although it is not the SU{Nf) x SU{Nf) of the continuum, it is a U{1). 
Thus there is only one Goldstone pion at flnite a, plus other non-degenerate 
pseudoscalar states whose mass goes to zero in the continuum limit (See Fig. ^ 
for an example of this.) They are preferred over Wilson fermions in situations 
in which the chiral properties of the fermions dominate the dynamics. They 
also cheaper to simulate than Wilson fermions, because there are less_variables. 
However, flavor and translational symmetry are all mixed togetherEj. 
(c) Chiral, undoubled, but not ultralocal 

These actions implement a modified version 113 of the chiral rotation 

5V = 75(l-^a^)V'; ^ i^il - ^aD)-f5 (41) 

which is sufficient to preserve all the interesting features of continuum chiral 
symmetry. An example of such an action is the "domain wall fermiontj." It is 
a variation on the idea that if you have a fermion coupled to a scalar field, and 
the scalar field interpolates between two minima (forms a soliton) , the fermion 
will develop a zero-energy chiral mode bound to the center of the soliton. Now 
we go into brane world, extend QCD into five dimensions, and put ourselves 
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and our four dimensional world on the kink. There is an anti-kink out there in 
the fifth dimension, and as long as it is far away the mode on the kink doesn't 
see the anti-kink and the 4-d theory on the kink is chiral. But if the anti-kink 
is too close (fifth dimension too small) the modes mix and chiral symmetry is 
broken. How close is "too close" is (yet) another engineering question. There 
are four dimensional analogs of this-think of integrating out the modes in 
the fifth dimension in favor of a tower of massive fermions, and get "overlap 
fermions til," or construct a low energy Wilsonian effective action from an 
underlying chiral theory and get "fixed point fermions t3." The bad feature 
is that these actions have couplings which reach out to many neighboring 
sites. Their strength drops exponentially with distance, so they are true local 
actions in the continuum limit, but they are very expensive to simulate. But 
stay tuned... 



3 Hadronic Matrix Elements from the Lattice 

One of the major goals of lattice calculations is to provide hadronic matrix 
elements which either test QCD or can be used as inputs to test the standard 
model. 



3.1 Generic Matrix Element Calculations 

Most of the matrix elements measured on the lattice are extracted from expec- 
tation values of local operators J{x) composed of quark and gluon fields. For 
example, if one wanted {0\J{x)\h) one could look at the two-point function 

Cjoit) =Y,{0\J{^,t)O{0,0m. (42) 
Inserting a complete set of correctly normalized momentum eigenstates 

and using translational invariance and going to large t gives 

CMt) = e-^^^M4Mm, (44) 
2mA 

A second calculation of 

Coo{t) = V(0|O(x,t)O(0,0)|0) ^ *MiJl^ (45) 
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is needed to extract (0| J|yl) by fitting two correlators with three parameters. 
Similarly, a matrix element {h\J\h') can be gotten from 

CAB{t,t')^Y.^Q\OA{t)J{x,t')OB{m)- (46) 

X 

by stretching the source and sink operators Oa and Ob far apart on the lattice, 
letting the lattice project out the lightest states, and then measuring and 
dividing out {Q\OA\h) and (0| Ob |^)- 

These lattice matrix elements are not yet the continuum matrix elements. 
Typically, one is interested in some matrix element defined with a particular 
regularization scheme. It is a generic feature of quantum field theory that an 
operator defined in one scheme {MS) will be a superposition of operators in 
another scheme (lattice). In principle, the superposition could be all possible 
operators. So generically an operator of dimension D will mix like 

The only restriction are symmetries: in a theory where parity is conserved a 
vector operator and an axial vector operator can't mix. This is relevant for 
lattice calculations because the symmetries of the lattice action are in general 
different from continuum symmetries. For example, the space-time symmetry 
of the lattice is given by the group of discrete rotations. A more serious source 
of mixing for light quark operators is the way lattice fermions treat chiral 
symmetry. Wilson-type fermions break chiral symmetry (even massless ones 
do so off-shell) and so nothing prevents mixing into "wrong chirality" operators. 

In Eq. |47| the "diagonal" term will contain the anomalous dimension of 
the continuum operator 

= l + -7-^(7rJog a^ + yl) + ... (48) 

(which cancels the mu-dependence of the coefficient function- C(/i)(/|0'^°"*|z)^ 
is independent of the renormalization point). In principle the leading log could 
be summed, but in practice we don't know how much of the constant term A 
should be absorbed into a change of scale of so they are just left there. 
The mixing terms to other dimension D operators die out in the continuum so 
they don't have any logs. There are also terms for mixing with higher dimen- 
sional operators, which give contributions proportional to positive powers of a. 
(These are usually benign.) One can also have mixing with lower dimensional 
operators, with contributions involving negative powers of a. (Four fermion 
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operators for Bk could mix with sd.) These are deadly. They must drop out 
in the continuum but it is a delicate business, since they look like they are 
growing as an inverse power of a. 

This is probably more than you wanted to know, but you need the Znm 's to 
produce numbers. People get them in a number of ways. Most straightforward 
is to compute them in perturbation theory, but lattice perturbation theory 
in terms of the bare coupling g(a) is not very convergent, and it is a long 
tricky story to do better. The culprit is the "tadpole graph." The lattice 
fermion-gauge field interaction is generically ip{x)U^{x)'ip{x + rhu) and U ~ 
1 + igaA^ — g^a^ /2A^ + . . .. The ^Aj^tp vertex, not present in any sensible 
continuum regularization, causes problems when the gluon forms a loop: the 
quadratic divergence from the loop integral combines with the to give a finite 
contribution-in fact, it is often the dominant contribution. In perturbation 
theory one must also choose the momentum scale in the .running coupling 
constant. There are reasonable choices for how to do thatBEa. 

Ofteii-pne can find Z„m's by forcing lattice observables to obey Ward 
identitiesE3. One can also play this game with quark propagators and vertices, 
by computing analogs of quatk vertices on the lattice and matching ones results 
to a continuum calculationcj. 

Besides, the Z's, there are other things that can go wrong. Most lattice 
actions break down when the quark mass gets heavy. The dispersion relation 
for Wilson or clover actions is E{p) = mi +p^/{2m2) and the quark magnetic 
moment is /i = l/wa with mi ^ TO2 ^ m^. The residue of the quark propa- 
gator at its pole is not 1/{2E) as in the continuum. What to do then is not 
obvious (meaning that lattice people fight over what to do). 

3.2 Heavy quark operators 

There are many lattice calculations of /s, /£), Bg^ and form factors for semilep- 
tonic decay. B — B mixing is parameterized by the ratio Xd — {AM),j^/r^^ 

= T,sp^ijQCDF{^)\V:,Vu?h{pi){\{B\hp{l - i,)dhp{l - l,)d\B)} 

(49) 

Experiment is on the left; theory on the right. Moving into the long equation 
from the left, we see many known (more or less) parameters from phase space 
integrals or perturbative QCD calculations, then a combination olCKM matrix 
elements, followed by a four quark hadronic matrix element E3. We would 
like to extract the CKM matrix element from the measurement of x^ (and 
its strange partner Xs). To do so we need to know the value of the object 
in the curly brackets, defined as 3/8Mm and parameterized as m%fg Bi,^ 
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where Bh^ is the so-caUed B-parameter, and fs is the B-meson decay constant, 
(0|67o75(i|i?) = fsiTT-B- Vacuum saturation suggests that Bb = 1- From the 
lattice one can try to get a real value. 

In Eq. ^ b{ij), the coefficient which runs the effective interaction down 
from the W-boson scale to the QCD scale /z, and the matrix element M{fi) 
both depend on the QCD scale. One often sees the renormalization group 
invariant quantities Mbd = i>{^)M},d(p) or Bm — b{p)BMi^J') quoted in the 
literature. 

Decay constants probe very simple properties of the wave function: in 
the nonrelativistic quark model /m = 4'{,Q)/^/itim, where ipi^) is the qq wave 
function at the origin. For a heavy quark (Q) light quark (g) system '0(0) 
should become independent of the heavy quark's mass as the Q mass goes to 
infinity, and in that limit one can show in QCD that ^tum^m approaches a 
constant. 

The decay constant is computed by combining a heavy quark and a light 
antiquark propagator into Eq. You might think it would be difficult to 
calculate directly on present day lattices with relativistic lattice fermions 
because the lattice spacing is much greater than the h quark's Compton wave- 
length (or the UV cutoff is below mh). But it is better to think of the lattice 
theory as an effective field theory for the low-momentum excitations in the pres- 
ence of additional high energy scales-the cutoff (inverse lattice spacing) and 
the heavy quark mass. As in any effective field theory, the effects of the short 
distance are lumped into coefficients of the effective theory E3. As a practical 
matter, one can use the good old clover action to do the calculations-it contains 
all the necessary operators. The bare mass has nothing direct to do with the 
results; one tunes it, monitoring the kinetic energy E(j)) = mi +p^/ (2m2) + . . ., 
and takes the hadron mass to be m2. 

Nonrelativistic QCD has also been discretizeiLand used to make very pre- 
cise calculations of the properties of quarkoniacJ. This formalism can also 
be used for the heavy quark (again as long as its momentum is small.) The 
"static" limit (infinite 6-quark mass) is often used as an additional point on 
the curve. Then one can try to extrapolate all the way from light quarks to 
heavies and get all the decay constants at once. 

I will show some pictures from the lattice decay constant of Ref. 0. These 
authors (my name is on it but I didn't do anything) did careful quenched sim- 
ulations at many values of the lattice spacing, which allows one to extrapolate 
to the continuum limit by brute force. They have also done a set of simula- 
tions which include light dynamical quarks, which should give some idea of the 
accuracy of the quenched approximation. 

Examples of results of Ref. l23 are shown in Figs. || and 0. The simulations 
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with dynamical fcrmions arc not as good quality as the quenched simulations: 
the lattice spacings are generally larger, the simulations all have two degenerate 
flavors (what about the strange quark), and the dynamical quark masses are 
still a bit large. We think that the Wilson results (crosses in Fig. ||(b) over es- 
timate the continuum result, and the clover action we are using underestimates 
it, but we also suspect that quenched /s is a bit too low. 



0.6 



Q used for-- 
used for- 
used for- 



CL=0.7£: 
CL=0.79 
CL=0.99 




0.0 



0.5 



1.0 



1/Mp (GeV ^) 



0. 



Figure 5: Pseudoscalar meson decay constant vs 1/M, from Refl 

Sonii has presented a summary of data from various collaborations, as 
of last winter. Again, there is a hint that the Nf = 2 results may be about 30 
MeV above the quenched results. 

Lattice calculations have been predicting quenched fn^ ~ 200 MeV for 
about twelve years. The central values have changed very little, while the 
uncertainties have decreased. Experimental determinations of fo^ all come in 
higher than the lattice results, though with large error bars. We need to do a 
good quality unquenched lattice calculation. 

Now back to the B parameter. On the lattice, one could measure the decay 
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> quenched Wilson 
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Figure 6: fs vs. a from Ref. 0, showing extrapolations to the continuum limit of quenched 
(a) and full (b) QCD data. The scale is set by fn = 132 MeV throughout. 



Table 1: Heavy-light decay constants and their ratios. 



Quantity 


Quenched {Nf = 0) 


Partially Unquenchcd {Nf = 2) 


/s/MeV 


170 ± 20 


200 ± 30 


/ss/McV 


190 ±20 


220 ± 30 


fo/MeV 


205 ± 20 


225 ± 30 




225 ± 20 


245 ± 30 


Ibs I Ib 


1.14 ± .06 


1.14 ± .06 


fos 1 fo 


1.10 ± .06 


1.10 ± .06 



constants and B parameter separately and combine them after extrapolation, 
or measure M directly. In principle the numbers should be the same, but 
in practice the first technique has produced better numbers so far. That is 
because the B parameter is measured as the ratio of a correlator with a four- 
fermion vertex to a product of two current-current correlators (see Fig. A 
lot of systematics cancel in the ratio. |— , |— , 

Many groups have visited this problem. Reviews by Draper and SoniEil 
quote a world summary, which I copy into Table ^. 

Semileptonic decays involve processes like Eqn. |4^. On the lattice, one 
just measures the matrix element of a current and fits it to the expected set 
of form factors-for B ■niv, for example, 

{^mV,W))=f+{q^)[p'+p- + (50) 

The best signals come when the momenta of the initial and final hadron are 
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Figure 7: Bk (shown) and Bb are computed rbv taking the ratio of four quark operator 
and two two-point functions. (Figure from Ref. e3.) 



Table 2: Summary of heavy-hght B-parameters. 





Quenched 


"Unqucnchcd" 


Bbs I Bbj 


.86(4)(8) 
1.00(1)(2) 
195 ± 25 MeV 

1.14 ± .06 


.86(4)(8) 
1.00(1)(2) 
230 ± 35 

1.14 ± .07 





small. Then the large B mass forces [q = lepton 4- momentum — pb — P-k) 
to be large. If the form factor is needed at ^ 0, a large extrapolation 
is needed, and there will be additional errors and model dependence in the 
answer. (Lattice people have no advantage over anyone else at guessing at 
functional forms.) However, finding Vub from experimental data only requires 
knowing the form factor at one value of q^. This should work so long as the 
experiment has enough data to measure the differential rate around that region 
of q^. Two recent approaches try to do this: UKQCD focussed on near the 
end-point or the zero-recoil region where the lattice data tends to be cleanest 
and heavy quark symmetry can be used. The FNAL group has measured 
B Div form factors at zero recoil L3. They have a clever technique from 
removing much of the lattice-to-continuum Z-factors by computing ratios of 
matrix elements, such as 

{D\cjob\B){B\hoc\D) , 

— — inwniT; , , m = \h+{v ■ v ^ 1)| (51) 

The denominators are just diagonal matrix elements of the charge density, 
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and they can easily be normalized. TheytJ are also computing semi-leptonic 
form factors for B — > iriiy and D Tr{K)£v, by concentrating directly on 
the differential decay spectrum in an interval with 0.4 < pT^/GeV < 0.8 thus 
avoiding the need for large extrapolation in g^. 



3.3 Kaon Matrix Elements 

Lattice calculations of kaon weak interaction matrix elements begin with the 
full Standard Model at high energies and use the operator product expansion, 
combined with the renormalization group, to construct a low-energy effective 
field theory valid at scales /i of a few GeV. The effective Hamiltonian basically 
reduces to a sum of four-fermion interactions 

10 

<^ = SE^*(^)0»(^) (52) 

^ i—O 

People have expended the most effort on, and have the best results for, Bk; 
there are some results on the A/ = 1/2 rule; and last, there is e'/e, with 
unreliable results so far. 
Bk 

The JLQCD collabO|Cation has the best results on Bk, from a calculation 
using staggered fcrmionfij. They have data from many lattice spacings and 
several choices for the lattice discretization of the operator. (See Fig. ^.) They 
find quenched Bk{MS, fi — 2 GeV) — 0.616(5). The main limitations of this 
result are quenching, plus the fact that the lattice calculations are actually done 
without SU{3) flavor breaking (the lattice "kaon" is a pseudoscalar made of 
degenerate quarks). These effects are believed&il to be 5-10 per cent corrections. 

JLQCD eII has also done a calculation with Wilson fermions. This was 
done not so much to get a number itself but to check the staggered result. The 
systematics are very different and the operator mixing is fierce due to the loss 
of chiral symmetry inherent in Wilson fermions. For example, S7;i(l — 75)^ ■ 
S7^(l — 75)d mixes with s^f^d ■ sj^d and that operator has a K — K matrix 
element ten times greater. 
A/= 1/2 Rule 

Can the lattice reproduce the experimentally observed factor of 22 between 
K'-' — ^ (7r7r)/=o and K'-^ {t^'^)i=2 amplitudes? The lattice calculations are 
difficult. In addition to graphs of Fig. |^, which are reasonably straightforward 
to compute, there are a host of other topologies, some of which involve com- 
puting propagators from many points on the lattice to many other points. But 
I think the reason there are so few lattice results is because all of the quantities 
of interest are scheme dependent and one must compute a lattice-to-continuum 
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Figure 8: Bk from staggered fermions, as a function of the-lattice spacing, for two different 
choices of lattice operators. (Figure and results from Ref. 



matching factor. In |addition, people don't calculate K —> tttt directly on the 
lattice; it is difhcultH to extract the phase shifts from the tttt final state in- 
teractions from lattice data (never mind trying to separate the two pions to 
|asymptotically great distances). Instead, they use chiral perturbation theory 
to relate K ^ tttt amplitudes to K n. In the case of the A/ = 3/2 
amplitude there is a factor of two change in the lattice result depending on 
whether tree level or one loop chiral perturbation theory is used. This is shown 
in Fig. |. 

— 

The only recent work I am aware of is by Pekurovsky and KilcupEj. The 
calculation is hard. The biggest operators Oq and Og in the nomenclature 
have opposite signs and nearly cancel. But the big problem is the scheme 
matching. In a perturbative calculation, we saw that (O)jj-g ~ Z{0)iatt and 
Z = 1 + ttsC -|- . . .. But in this equation, what scheme is used to compute 
as , and what is the scale q* at which as is evaluated? Pekurovsky and Kilcup 
found that their numbers for one operator, Og, shifted by a factor of 4 when 
they were converted into MS, and the factor of 4 could become a factor of 30 
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Figure 9: A/ = 3/2 K — > -kit amplitude with (fancy symbols) and without (plain symbols) 
one-loon— corrections of quenched chiral perfr»rbation theory. Data are crosses and fancy 
crosses, c3; diamonds and fancy diamonds, CJ- iBjata are plotted as a function of lattice 
meson mass. This figure (and much else) is from r\ 



(or worse) as q* was varied from it /a to 1/a. They attempt to guesstimate 
numbers but since they say plainly that they should be used "with extreme 
caution" I won't quote them. A nonperturbative approach to matching is 
clearly needed but does not exist yet. 



4 Conclusions 

What about the future? Matrix elements are at the end of a long chain involv- 
ing a large set of both simulation and physics issues. They are the most com- 
plicated corner of the lattice game. If all you want are the numbers, Moore's 
law says that computer speed doubles every eighteen months, and statistics is 
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"v/TV, so error bars will fall by a factor of 2 every three years for everything 
we know how to do today and will learn nothing new about how to do better 
in the future. And there are many projects and proposals to build clusters 
or dedicated supercomputers at a cost which is "chicken feed" compared to 
the experimental program. This will enable us to begin to chip away at the 
biggest systematic in all the calculations I have shown here-the neglect of the 
quenched approximation. 

But new hardware is not really where the action is. It is merely "enabling 
technology," so we can make mistakes faster, learn more about the physics, 
and test new ideas. 

The main bottleneck to progress on hadronic matrix elements is just that 
these calculations are complex and have many parts. Some of us (me, let's not 
be shy) think that better discretization algorithms will help. The problem with 
that approach is that many pieces of the puzzle have to be determined from 
scratch: learning how to optimize the new algorithm, testing spectroscopy, 
computing the Z's. This takes a couple of years, if the inventor of the algorithm 
doesn't get tired first. Others of us prefer to live with poorer algorithms (which 
have already been well calibrated) and try to tweak the parts of them which 
work the worst. The simulations still take a couple of years. Believe it or not, 
even though lattice QCD is a mature field, there are still many questions about 
QCD which lattice people do not know how to answer, and an outsider might. 
Maybe you would enjoy thinking about them. 
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